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Abstract — Parallel least mean square-partial parallel interfer- 
ence cancelation (PLMS-PPIC) is a partial interference cance- 
lation which employs adaptive multistage structure [1]. In this 
algorithm the channel phases for all users are assumed to be 
known. Having only their quarters in (0, 2it), a modified version 
of PLMS-PPIC is proposed in this paper to simultaneously esti- 
mate the channel phases and the cancelation weights. Simulation 
examples are given in the cases of balanced, unbalanced and 
time varying channels to show the performance of the modified 
PLMS-PPIC method. 

I. Introduction 

The multiple access interferences (MAI) is the root of user 
limitation in CDMA systems [2], [5]. The parallel least mean 
square-partial parallel interference cancelation (PLMS-PPIC) 
method is a multiuser detector for code division multiple 
access (CDMA) receivers which reduces the effect of MAI in 
bit detection. In this method and similar to its former versions 
like LMS-PPIC [6] (see also [7]), a weighted value of the 
MAI of other users is subtracted before making the decision 
for a specific user in different stages [1]. In both of these 
methods, the normalized least mean square (NLMS) algorithm 
is engaged [8]. The m"^ element of the weight vector in each 
stage is the true transmitted binary value of the m"^ user 
divided by its hard estimate value from the previous stage. 
The magnitude of all weight elements in all stages are equal 
to unity. Unlike the LMS-PPIC, the PLMS-PPIC method tries 
to keep this property in each iteration by using a set of NLMS 
algorithms with different step-sizes instead of one NLMS 
algorithm used in LMS-PPIC. In each iteration, the parameter 
estimate of the NLMS algorithm is chosen whose element 
magnitudes of cancelation weight estimate have the best match 
with unity. In PLMS-PPIC implementation it is assumed that 
the receiver knows the phases of all user channels. However in 
practice, these phases are not known and should be estimated. 
In this paper we improve the PLMS-PPIC procedure [1] in 
such a way that when there is only a partial information of the 
channel phases, this modified version simultaneously estimates 
the phases and the cancelation weights. The partial information 
is the quarter of each channel phase in (0, 2tt). 

The rest of the paper is organized as follows: In section HJ 
the modified version of PLMS-PPIC with capability of channel 
phase estimation is introduced. In section |III] some simulation 



examples illustrate the results of the proposed method. Finally 
the paper is concluded in section |lVl 

II. Multistage Parallel Interference 
Cancelation: Modified PLMS-PPIC Method 

We assume M users synchronously send their symbols 
ai, a2, • • • , anj via a base-band CDMA transmission system 
where G {—1, 1}. The m*'' user has its own code Pm{-) 
of length TV, where Pm{n) G { — 1,1}, for all ri. It means 
that for each symbol N bits are transmitted by each user and 
the processing gain is equal to A^. At the receiver we assume 
that perfect power control scheme is applied. Without loss of 
generality, we also assume that the power gains of all channels 
are equal to unity and users' channels do not change during 
each symbol transmission (it can change from one symbol 
transmission to the next one) and the channel phase (j>„i of 
m*'' user is unknown for all m ~ 1,2, ■ ■ ■ , M (see [1] for 
coherent transmission). According to the above assumptions 
the received signal is 



M 



r{n) 



p]<t>r, 



'Pm{n) + v{n), n = l,2,---,N, (1) 



where v{n) is the additive white Gaussian noise with zero 
mean and variance . Multistage parallel interference can- 
celation method uses a\^^ , a'^^'^ , ■ ■ ■ , a^J^, the bit estimates 
outputs of the previous stage, s— 1, to estimate the related MAI 
of each user. It then subtracts it from the received signal r{n) 
and makes a new decision on each user variable individually 



to make a new variable set a 



a%i for the current 



stage s. Usually the variable set of the first stage (stage 0) is 
the output of a conventional detector The output of the last 
stage is considered as the final estimate of transmitted bits. In 
the following we explain the structure of a modified version 
of the PLMS-PIC method [1] with simultaneous capability of 
estimating the cancelation weights and the channel phases. 
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1. Define 



Assume a\n. £ { — 1, 1} is a given estimate of from 



,,(^-1) 



(2) 



From ([T]) and (|2]i we have 



r{n) = ^ w'^^ai^ ^^Pm{n) + v{n). 



(3) 



Define 



W = K,«^2,- (4a) 
X%n) = [al^-SiN,4'"S2N,---,4rSMN]?4b) 

where T stands for transposition. From equations (l4al i and 
J, we have 



r{n) = W 



(5) 



Given the observations {r{n), X*(n)}^^]^, in modified PLMS- 
PPIC, like the PLMS-PPIC [1], a set of NLMS adaptive 
algorithm are used to compute 



W^{N) = [wl{N),w^,{N),---,wlAN)f, 



(6) 



which is an estimate of W after iteration N. To do so, from 
we have 



which is equivalent to 



M 



11 =0. 



(7) 



(8) 



We divide 4* = |^0, 1 — y %/ , a sharp range for /i (the step- 
size of the NLMS algorithm) given in [9], into L subintervals 
and consider L individual step-sizes 8 = {/ii, /i2, ■ • • 7 Ml}' 

where = — , /i2 = 2/ii, • • •, and /i^ — L^i. In each 
stage, L individual NLMS algorithms are executed (fii is the 
step-size of the algorithm). In stage s and at iteration n, 
if W^{n) = [wf fc, • • • , w]^ 1.]'^, the parameter estimate of the 
k^^ algorithm, minimizes our criteria, then it is considered as 
the parameter estimate at time iteration n. In other words if 
the next equation holds 



arg mm 




(9) 



where Wf{n) = W'{n - 1) + Mi j, y^ffp e(n), / = 
1,2,- l,i and = {Wf{n),---,Wl{n)}, 
then we have ^{n) = W^{n), and therefore all other 
algorithms replace their weight estimate by W^{n). At time 
instant n = N, this procedure gives ^{N), the final estimate 
of W'\ as the true parameter of stage s. 

Now consider R = (0, 2tt) and divide it into four equal parts 
i?i = (0, f ), i?2 = (|,^), R3 = (tt, \) and i?4 = (f ,27r). 
The partial information of channel phases (given by the 
receiver) is in a way that it shows each 0„i (m = 1, 2, • • • , M) 
belongs to which one of the four quarters i?;, i ~ 1,2,3,4. 
Assume = [wl{N),wl{N), - ■ ■ ,wl^{N)]^ is the 



weight estimate of the modified algorithm PLMS-PPIC at time 
instant N of the stage s. From equation (|2]i we have 

(s-l) 



4>m = ^i- 



We estimate 



by 



where 

(s-l) 



-win)- 



Because — — = 1 or —1, we have 



Z<(7V) 
±^ + Z<(iV) if 



-1 



(10) 



(11) 



(12) 



Hence 0f„ G P« = {^<(iV), ^<(iV) + tt, Z<(^) - tt}. 
If wf^{N) sufficiently converges to its true value wf^^, the 
same region for and (pm is expected. In this case only 
one of the three members of P'' has the same region as 
For example if 0m S (0, §), then 0^ G (0, -j) and therefore 
only ZWjj^{N) or Zw!^{N) + tt or Zw!l^{N) — n belongs 
to (0, ^). If, for example, Zwfj^{N) + tt is such a member 
between all three members of P*, it is the best candidate for 
phase estimation. In other words, 

<^™«C = Z<(Af)+7r. 

We admit that when there is a member of in the quarter 
of (pm, then u)^(iV) converges. What would happen when 
non of the members of P** has the same quarter as 0m? This 
situation will happen when the absolute difference between 
Zwf^{N) and (pm is greater than tt. It means that wf^{N) 
has not converged yet. In this case where we can not count 
on w!l^{N), the expected value is the optimum choice for the 
channel phase estimation, e.g. if (pm G (Oi f ) then j is the 
estimation of the channel phase (/),„, or if (/),„ S (f j'"') then 
^ is the estimation of the channel phase (pm- The results of 
the above discussion are summarized in the next equation 

Z<(7V) if Z<(7V),0m eP,, z = 1,2,3,4 

Zw^m{N)+TT if Zu)^(iV) +7r,0,„ G P„ i = 1,2,3,4 

Zw'^nW-n if Zw'^{N)~^T,(j)m&R^, i = 1,2,3,4 

il^lh±l^ if PmeR^, Zwf,,{N),Zwf^{N)±niR„ 

Having an estimation of the channel phases, the rest of the 
proposed method is given by estimating as follows: 




(13) 



{N)c^t'"Pr,An)- (14) 



1 — l,m 



The inputs of the first stage {oi^i}m=i (needed for computing 
X^{n)) are given by 



N 



am = sign < real I ^ ?'('^)e ^'^"^Pm{n) 



(15) 



TABLE I 

Channel phase estimate of the first user (example[TJ 



TABLE II 

Channel phase estimate of the first user (example^ 
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Assuming ipm G Ri, then 



Table HU] shows the structure of the modified PLMS-PPIC 
method. It is to be notified that 

• Equation ( fTSl l shows the conventional bit detection 
method when the receiver only knows the quarter of 
channel phase in (0, 27r). 

• With L = 1 (i.e. only one NLMS algorithm), the modified 
PLMS-PPIC can be thought as a modified version of the 
LMS-PPIC method. 

In the following section some examples are given to illus- 
trate the effectiveness of the proposed method. 

III. Simulations 

In this section we have considered some simulation exam- 
ples. Examples [T][3] compare the conventional, the modified 
LMS-PPIC and the modified PLMS-PPIC methods in three 
cases: balanced channels, unbalanced channels and time vary- 
ing channels. In all examples, the receivers have only the 
quarter of each channel phase. Example [T]is given to compare 
the modified LMS-PPIC and the PLMS-PPIC in the case of 
balanced channels. 

Example 1: Balanced channels: Consider the system 
model (O in which M users synchronously send their bits 
to the receiver through their channels. It is assumed that each 
user's information consists of codes of length N . It is also as- 
sumd that the signal to noise ratio (SNR) is OdB. In this exam- 
ple there is no power-unbalanced or channel loss is assumed. 
The step-size of the NLMS algorithm in modified LMS-PPIC 

method is = 0.1(1 — \J~^^^) and the set of step-sizes 
of the parallel NLMS algorithms in modified P LMS- PPIC 
method are 6 = {0.01, 0.05, 0.1, 0.2, • • • , 1}(1 - \J i.e. 

Ati = 0.01(1-^^), • • • , = 0.2(l-y^), • • • ,//i2 = 

(1 - \f^^)- Figure [T] illustrates the bit error rate (BER) 
for the case of two stages and for = 64 and N — 256. 
Simulations also show that there is no remarkable difference 
between results in two stage and three stage scenarios. Table U 
compares the average channel phase estimate of the first user 
in each stage and over 10 runs of modified LMS-PPIC and 
PLMS-PPIC, when the the number of users is M = 15. 

Although LMS-PPIC and PLMS-PPIC, as well as their 
modified versions, are structured based on the assumption of 
no near-far problem (examples |2] and O, these methods and 



especially the second one have remarkable performance in the 
cases of unbalanced and/or time varying channels. 

Example 2: Unbalanced channels: Consider example [T| 
with power unbalanced and/or channel loss in transmission 
system, i.e. the true model at stage s is 

AI 

r{n) = Pmw'^al^l'^^ Cmin) + v{n), (17) 

m— 1 

where < /3™ < 1 for all 1 < m < M. Both the LMS- 
PPIC and the PLMS-PPIC methods assume the model Q, and 
their estimations are based on observations {r{n), X'^ (n)}, 
instead of {r{n),GX^{n)}, where the channel gain matrix 
is G = diag(/3i, /32, • • • , (3m)- In this case we repeat example 
[T] We randomly get each element of G from [0, 0.3]. Figure |2] 
illustrates the BER versus the number of users. Table |ll] 
compares the channel phase estimate of the first user in each 
stage and over 10 runs of modified LMS-PPIC and modified 
PLMS-PPIC for M = 15. 

Example 3: Time varying channels: Consider example [T] 
with time varying Rayleigh fading channels. In this case 
we assume the maximum Doppler shift of 40HZ, the 
three-tap frequency-selective channel with delay vector of 
{2 X 10-^2.5 X 10-^3 X lO-'^jsec and gain vector of 
{— 5, — 3, — lOjdB. Figure [3] shows the average BER over all 
users versus M and using two stages. 

IV. Conclusion 

In this paper, parallel interference cancelation using adaptive 
multistage structure and employing a set of NLMS algorithms 
with different step-sizes is proposed, when just the quarter 
of the channel phase of each user is known. In fact, the 
algorithm has been proposed for coherent transmission with 
full information on channel phases in [1]. This paper is a 
modification on the previously proposed algorithm. Simu- 
lation results show that the new method has a remarkable 
performance for different scenarios including Rayleigh fading 
channels even if the channel is unbalanced. 
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Fig. 1. The BER of the conventional, the modified LMS-PPIC and the 
modified PLMS-PPIC methods versus the system load in balanced channel, 
using two stages for N = 64 and N = 256. 
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Fig. 2. The BER of the conventional, the modified LMS-PPIC and the 
modified PLMS-PPIC methods versus the system load in unbalanced channel, 
using two stages for N = 6A and N = 256. 



Fig. 3. The BER of the conventional, the modified LMS-PPIC and 
the modified PLMS-PPIC methods versus the system load in time varying 
Rayleigh fading channel, using two stages for A'^ = 64 and A'^ = 256. 
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TABLE III 

The procedure of the modified PLMS-PPIC method 



for m = 1,2,- • •,A'/ 
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i^m 4 



TV 
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for s = 1, 2, • • • , 5 



o 



Oh 



Oh 



for n = 1,2, •••,iV 



for m = 1, 2, • 



.M 



for 771 = 1, 2, • • • , M 



e(7i) = r(77) - M/"'^(77 - l)X''(r7) 
Z{n) = ||j^I(g{. e(77) 
min = 00, 1 = 1 
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W'{n) = Wfin) 
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if E Ifc.feWI -1| <min: 



m— 1 
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/ = k 
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